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Abstract 

The Higgs boson search has shifted from LEP2 to the Tevatron and wih subse- 
Q^' quently move to the LHC. The current hmits from the Tevatron and the prospective 

^ • sensitivities at the LHC are often interpreted in specific MSSM scenarios. For heavy 

D , Higgs boson production and subsequent decay into bb or r+r", the present Tevatron 

' data allow to set limits in the M_4-tan /3 plane for small Ma and large tan /3 values. 

Similar channels have been explored for the LHC, where the discovery reach extends 
^ . to higher values of and smaller tan (3. Searches for MSSM charged Higgs bosons, 

^ I produced in top decays or in association with top quarks, have also been investigated 

at the Tevatron and the LHC. We analyze the current Tevatron limits and prospective 
LHC sensitivities. We discuss how robust they are with respect to variations of the 
other MSSM parameters and possible improvements of the theoretical predictions for 
Higgs boson production and decay. It is shown that the inclusion of super symmetric 
radiative corrections to the production cross sections and decay widths leads to impor- 
tant modifications of the present limits on the MSSM parameter space. The impact 
on the region where only the lightest MSSM Higgs boson can be detected at the LHC 
is also analyzed. We propose to extend the existing benchmark scenarios by including 
additional values of the higgsino mass parameter /i. This affects only slightly the search 
channels for a SM-like Higgs boson, while having a major impact on the searches for 
non-standard MSSM Higgs bosons. 
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1 Introduction 



Disentangling the mechanism that controls electroweak symmetry breaking is one of the main 
tasks of the current and next generation of colliders. Among the most studied candidates 
in the literature are the Higgs mechanism within the Standard Model (SM) or within the 
Minimal Supcrsymmctric Standard Model (MSSM). Contrary to the SM, two Higgs doublets 
are required in the MSSM, resulting in five physical Higgs boson degrees of freedom. In the 
absence of explicit CP- violation in the soft supersymmetry-breaking terms these are the light 
and heavy CV-cvcn Higgs bosons, h and H, the CP-odd Higgs boson. A, and the charged 
Higgs boson, H^. The Higgs sector of the MSSM can be specified at lowest order in terms 
of Mz, Ma, and tan/3 = V2/V1, the ratio of the two Higgs vacuum expectation values. The 
masses of the CP-even neutral Higgs bosons and the charged Higgs boson can be calculated, 
including higher-order corrections, in terms of the other MSSM parameters. 

After the termination of LEP in the year 2000 (the close-to-final LEP results can be found 
in Refs. [1,2]), the Higgs boson search has shifted to the Tevatron and will later be continued 
at the LHC. Due to the large number of free parameters, a complete scan of the MSSM 
parameter space is too involved. Therefore the search results at LEP have been interpreted [2] 
in several benchmark scenarios [3,4]. Current analyses at the Tevatron and investigations of 
the LHC potential also have been performed in the scenarios proposed in Refs. [3,4]. The 
m™'^'^ scenario has been used to obtain conservative bounds on tan j3 for fixed values of the 
top-quark mass and the scale of the supcrsymmctric particles [5]. Besides the m™^^ scenario 
and the no-mixing scenario, where a vanishing mixing in the stop sector is assumed, the 
"small cteff" scenario and the "gluophobic Higgs scenario" have been investigated [6] . While 
the latter one exhibits a strong suppression of the ggh couphng over large parts of the Ma- 
tan f3 parameter space, the small ctcff scenario has strongly reduced couplings of the light 
CP-even Higgs boson to bottom-type fermions up to Ma ^ 350 GeV. These scenarios are 
conceived to study particular cases of challenging and interesting phenomenology in the 
searches for the SM-hke Higgs boson, i.e. mostly the hght CP-even Higgs boson. 

The current searches at the Tevatron are not yet sensitive to a SM-like Higgs in the 
mass region allowed by the LEP exclusion bounds [1,2]. On the other hand, scenarios with 
enhanced Higgs boson production cross sections can be probed already with the currently 
accumulated luminosity. Enhanced production cross sections can occur in particular for low 
Ma in combination with large tan /3 due to the enhanced couplings of the Higgs bosons to 
down-type fermions. The corresponding limits on the Higgs production cross section times 
branching ratio of the Higgs decay into down-type fermions can be interpreted in MSSM 
benchmark scenarios. Limits from Run II of the Tevatron have recently been published for 
the following channels [7-9] (here and in the following denotes all three neutral MSSM 
Higgs bosons, (f) — h, H, A): 



The obtained cross section limits have been interpreted in the m™'^'' and the no-mixing 
scenario with a value for the higgsino mass parameter oi /i — —200 GeV [7] and /j, — 



bh(j), (j) ^ hh (with one additional tagged h jet), 
pp ^ (f) ^ T~^T~ (inclusive), 
pp^tt^ H^W"^ bb, . 



(1) 
(2) 
(3) 
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±200 GeV [8]. In these scenarios for ~ 100 GcV the limits on tan/5 are tan/9 < 50. 

In this article, we investigate the dependence of the CDF and DO exclusion bounds in 
the M^-tan (3 plane on the parameters entering through the most relevant super symmetric 
radiative corrections in the theoretical predictions for Higgs boson production and decay 
processes. We will show that the bounds obtained from the bb(f), (p ^ bb channel depend 
very sensitively on the radiative corrections affecting the relation between the bottom quark 
mass and the bottom Yukawa coupling.^ In the channels with t~^t~ final states, on the other 
hand, compensations between large corrections in the Higgs production and the Higgs decay 
occur. In this context we investigate the impact of a large radiative correction in the gg ^ (f) 
production process that had previously been omitted. 

In order to reflect the impact of the corrections to the bottom Yukawa coupling on the 
exclusion bounds we suggest to supplement the existing m^^^ and no-mixing scenarios, 
mostly designed to search for the hght CP-even MSSM Higgs boson, h, with additional 
values for the higgsino mass parameter In fact, varying the value and sign of /i, while 
keeping fixed the values of the ghiino mass and the common third generation squark mass 
parameter Msusy, demonstrates the effect of the radiative corrections on the production 
and decay processes. The scenarios discussed here are designed specifically to study the 
MSSM Higgs sector without assuming any particular soft supersymmetry-breaking scenario 
and taking into account constraints only from the Higgs boson sector itself. In particular, 
constraints from requiring the correct cold dark matter density, BR(6 sj) or {g — 2)^, 
which depend on other parameters of the theory, are not crucial in defining the Higgs boson 
sector, and may be avoided. However, we also include a brief discussion of the "constraincd- 
^max)5 scenario, which in the case of minimal flavor violation and positive values of fi leads 
to a better agreement with the constraints from BR(6 — > 57). 

We also study the non-standard MSSM Higgs boson search sensitivity at the LHC, fo- 
cusing on the processes pp H/A + X, H/A — r+r" and pp tH^ + X, tu^, 
and stress the relevance of the proper inclusion of supersymmctric radiative corrections to 
the production cross sections and decay widths. We show the impact of these corrections by 
investigating the variation of the Higgs boson discovery reach in the benchmark scenarios for 
different values of /x. In particular, we discuss the resulting modification of the parameter 
region in which only the light CP-even MSSM Higgs boson can be detected at the LHC. 

The paper is organized as follows: Section 2 gives a summary of the most relevant 
supersymmctric radiative corrections to the Higgs boson production cross section and decay 
widths, while also introducing our notation. In section 3 we analyze the impact of these 
radiative corrections on the current Tevatron limits in the large tan f3 region, as well as on 
the future LHC reach for the heavy, non-standard, MSSM Higgs bosons. In section 4, based 
on the results of section 3, we propose an extension of the existing benchmark scenarios. 
The conclusions are presented in section 5. 

^We concentrate here on the effects of supersymmetric radiative corrections. For a recent account of 
uncertainties related to parton distribution functions, see e.g. Ref. [10]. 
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2 Predictions for Higgs boson production and decay 
processes 

2.1 Notation and renormalization 

The tree-level values for the CP-even Higgs bosons of the MSSM, rrih and mn, are determined 
by tan/3, the CP-odd Higgs-boson mass Ma, and the Z boson mass Mz- The mass of the 
charged Higgs boson, Mh±, is given in terms of Ma and the W boson mass, Myy. Beyond 
the tree-level, the main correction to the Higgs boson masses stems from the t/t sector, and 
for large values of tan/5 also from the b/b sector. 

In order to fix our notations, we list the conventions for the inputs from the scalar top and 
scalar bottom sector of the MSSM: the mass matrices in the basis of the current eigenstates 
ii, Ir and 6^, Ir are given by 

^2 _ ( Ml+ml + co,2H\-lsl)Ml m,X, \ 

~ 1^ mtXt M^+m2 + |cos2/?4M| j' 

■^^5 ~ 1^ m.X, Ml^ + ml-\cos2l3slMl j' 

where 

rritXt = mt{At — ^ cot (3 ) , rribXf, = rrih {Af, — fitanP ). (6) 

Here At denotes the trilinear Higgs-stop coupling. Ah denotes the Higgs-sbottom coupling, 
and II is the higgsino mass parameter. 

SU(2) gauge invariance leads to the relation 

Mi,=M-,^. (7) 

For the numerical evaluation, a convenient choice is 

Mt^ = M^^ =Mt^ = M^^ =: MsusY. (8) 

Concerning analyses for the case where M^^ ^ Mi^ ^ ^Iri ^-S- R-efs. [11,12]. It has 
been shown that the upper bound on rrih obtained using eq. (jH)) is the same as for the more 
general case, provided that Msusy is identified with the heaviest mass of M^^, M^^, M^^ [12]. 
We furthermore use the short-hand notation 

Ml:=Ml^^^ + ml . (9) 

Accordingly, the most important parameters for the corrections in the Higgs sector are 
mt, Msusy, -^t, and X^. The Higgs sector observables furthermore depend on the SU(2) 
gaugino mass parameter, M2. The other gaugino mass parameter. Mi, is usually fixed via 
the GUT relation ^ 

Ml = |%M2. (10) 

At the two-loop level also the gluino mass, m^, enters the predictions for the Higgs-boson 
masses. 
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Corrections to the MSSM Higgs boson sector have been evaluated in several approaches. 
The status of the available corrections to the masses and mixing angles in the MSSM Higgs 
sector (with real parameters) can be summarized as follows. For the one-loop part, the 
complete result within the MSSM is known [13-16]. The by far dominant one-loop contri- 
bution is the 0{at) term due to top and stop loops {at = /ij/(47r), ht being the top-quark 
Yukawa coupling). Concerning the two- loop effects, their computation is quite advanced 
and has now reached a stage such that all the presumably dominant contributions are 
known [12, 17-31]. They include (evaluated for vanishing external momenta) the strong 
corrections, usually indicated as 0{atas), and Yukawa corrections, 0{a^), to the dominant 
one-loop 0{at) term, as well as the strong corrections to the bottom/ sbottom one-loop 0{ab) 
term (aj, = /i^/(47r)), i.e. the 0{ai,as) contribution. The latter can be relevant for large val- 
ues of tan j3. For the (s)bottom corrections the all-order resummation of the tan /9-enhanced 
terms, 0{ab{as tan f3)'^), has also been computed. Recently the 0{atab) and 0{al) cor- 
rections have been obtained. The remaining theoretical uncertainty on the light CP-even 
Higgs boson mass has been estimated to be below ~ 3 GeV [32,33]. The above calculations 
have been implemented into public codes. The program FeynHiggs [34-37] is based on the 
results obtained in the Feynman-diagrammatic (FD) approach [12,21,31,32]. It includes all 
the above corrections. The code CPsuperH [38] is based on the renormalization group (RG) 
improved effective potential approach [18, 19,39]. For the MSSM with real parameters the 
two codes can differ by up to ~ 4 GeV for the light CP-even Higgs boson mass, mostly due 
to formally subleading two-loop corrections that are included only in FeynHiggs. For the 
MSSM with complex parameters the phase dependence at the two-loop level is included in a 
more advanced way [40] in CPsuperH, but, on the other hand, CPsuperH does not contain all 
the subleading one-loop contributions that are included [41,42] in FeynHiggs. Most recently 
a full two-loop effective potential calculation (including even the momentum dependence for 
the leading pieces) has been published [43]. However, no computer code is publicly available. 
In the following we will concentrate on the MSSM with real parameters. 

It should be noted in this context that the FD result has been obtained in the on-shell 
(OS) renormalization scheme, whereas the RG result has been calculated using the MS 
scheme; see Refs. [39,44] for a detailed comparison. Owing to the different schemes used in 
the FD and the RG approach for the renormahzation in the scalar top sector, the parameters 

Xt and MsusY arc also scheme-dependent in the two approaches. This difference between 
the corresponding parameters has to be taken into account when comparing the results of 
the two approaches. In a simple approximation the relation between the parameters in the 
different schemes is at 0{as) given by [39] 



where in the terms proportional to it is not necessary to distinguish between MS and on- 
shell quantities, since the difference is of higher order. The MS top-quark mass, m^^(mt) = 
rfit, is related to the top-quark pole mass, mf^ = rrit, in 0{as) by 



Ms' 



■2,MS 



3 TT 



(11) 




(12) 



mt = 



1+ 3^a,(mt) 



(13) 
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While the resulting shift in the parameter Msusy turns out to be relatively small in gen- 
eral, sizable differences can occur between the numerical values of Xt in the two schemes, 
see Refs. [12,39]. For this reason we specify below different values for Xi within the two 
approaches. 



2.2 Leading effects from the bottom/sbottom sector 

The relation between the bottom-quark mass and the Yukawa coupling /i^, which con- 
trols also the interaction between the Higgs fields and the sbottom quarks, reads at lowest 
order = hi^vi. This relation is affected at one-loop order by large radiative correc- 
tions [28-30,45], proportional to h}yV2i in general giving rise to tan /5-enhanced contributions. 
These terms proportional to ^2, often called threshold corrections to the bottom mass, are 
generated either by gluino-sbottom one-loop diagrams (resulting in 0{ahas) corrections to 
the Higgs masses), or by chargino-stop loops (giving 0{ai,at) corrections). Because the 
tan /3-enhanced contributions can be numerically relevant, an accurate determination of hb 
from the experimental value of the bottom mass requires a resummation of such effects to 
all orders in the perturbative expansion, as described in Refs. [29,30]. 

The leading effects are included in the effective Lagrangian formalism developed in 
Ref. [29]. Numerically this is by far the dominant part of the contributions from the sbottom 
sector (see also Refs. [26,27,31]). The dominant contributions arise from the loop-induced 
coupling of Hu (the Higgs field that couples at the tree-level to up-type fermions only) to 
the down-type fermions. The effective Lagrangian is given by 



^ 9 rrib 



2Mw 1 + Afc 



tan P Ai b-f^b + V2 Vtb tan j3 H^tLb 



R 



/ sm a , cos a \ - , 

+ --A,——]hbLbR 

\ cos p smp J 

/cosa . sina \ , 1 , 

- ^ + Afc— T HbLbR + h.c. . 14 

ycosp smp J J 

Here nTf, denotes the running bottom quark mass including SM QCD corrections. In the 
numerical evaluations obtained with FeynHiggs below we choose mb = m{,(mt) ^ 2.97 GeV. 
The prefactor 1/(1 -|- A^) in eq. ()14|1 arises from the resummation of the leading corrections 
to all orders. The additional terms ~ A;, in the hbb and Hbb couplings arise from the mixing 
and coupling of the "other" Higgs boson, H and h, respectively, to the b quarks. 

As explained above, the function A;, consists of two main contributions, an 0{as) cor- 
rection from a sbottom-gluino loop and an 0{at) correction from a stop-higgsino loop. The 
explicit form of A^ in the limit of Ms ^ rrit and tan (3^1 reads [28] 

Ab = ^m^/i tan/? x I (mi ,mr,mg) + ^ At tan P x I{mi mi^) . (15) 
The function / is given by 

1 



2 52 ^2^ 
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The large b — g loops are resummed to all orders of (a^ tan /?)"■ via the inclusion of A;, [28-30]. 
The leading electroweak contributions are taken into account via the second term in eq. ()15|) . 

For large values of tan /3 and the ratios of fimg/M^^^y and fiAt/M^^^y, the correction 
can become very important. Considering positve values of At and rrig, the sign of the A;, term 
is governed by the sign of /i. Cancellations can occur if At and rrig have opposite signs. For 
fj,,mg,At > the A;, correction is positive, leading to a suppression of the bottom Yukawa 
coupling. On the other hand, for negative values of A;,, the bottom Yukawa coupling may 
be strongly enhanced and can even acquire non-perturbative values when A;, — > —1. 

2.3 Impact on Higgs production and decay at large tan/5 

Higgs-boson production and decay processes at the Tevatron and the LHC can be affected 
by different kinds of large radiative corrections. The SM and MSSM corrections to the 
production channel gg ^ (p have been calculated in Refs. [46,47], SM corrections to the bb(j) 
channel have been evaluated in Refs. [48-50]. Higgs decays to bb and to t~^t~ within the 
SM and MSSM have been evaluated including higher-order corrections in Refs. [16,42,51]. 
Besides the process-specific corrections to the production and decay processes, large Higgs- 
boson propagator contributions have an impact on the Higgs-boson couplings. For large 
tan/3 the supersymmetric radiative corrections to the bottom Yukawa coupling described 
above become particularly important [52]. Their main effect on the Higgs-boson production 
and decay processes can be understood from the way the leading contribution A;, enters. 
In the following we present simple analytic approximation formulae for the most relevant 
Higgs-boson production and decay processes. They are meant for illustration only so that the 
impact of the Af, corrections can easily be traced (for a discussion of possible enhancement 
factors for MSSM Higgs-boson production processes at the Tevatron and the LHC, see also 
Refs. [53,54]). In our numerical analysis below, we use the full result from FeynHiggs rather 
than the simple formulae presented in this section. No relevant modification to these results 
would be obtained using CPsuperH. 

We begin with a simple approximate formula that represents well the MSSM parametric 
variation of the decay rate of the CP-odd Higgs boson in the large tan/3 regime. One 
should recall, for that purpose, that in this regime the CP-odd Higgs boson decays mainly 
into r-leptons and bottom-quarks, and that the partial decay widths are proportional to 
the square of the Yukawa couplings evaluated at an energy scale of about the Higgs boson 
mass. Moreover, for Higgs boson masses of the order of 100 GeV, the approximate relations 
mb(M^)^ ~ 9 GeV^, and mr(MA)^ — 3 GeV^ hold. Hence, since the number of colors is 
Nc = 3, for heavy supersymmetric particles, with masses far above the Higgs boson mass 
scale, one has 

BR{A bb) 
BR{A -> T+T-) 

On the other hand, the production cross section for a CP-odd Higgs boson produced 
in association with a pair of bottom quarks is proportional to the square of the bottom 



2 > (17) 

(1 + A,)V9' 

^1±^ . (18) 
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Yukawa coupling and therefore is proportional to tan^ /?/(! + A^)^. Also in the gluon fusion 
channel, the dominant contribution in the large tan/3 regime is governed by the bottom 
quark loops, and therefore is also proportional to the square of the bottom Yukawa coupling. 
Since the top-quark coupling is suppressed by either loop-corrections or inverse powers of 
tan /3, the leading top-quark correction arises from interference terms between the top-quark 
and bottom-quark loop diagrams. We have checked that these interference terms lead to 
corrections smaller than one percent (a few percent) for values of tan/3 ~ 50 (20). These 
corrections are small, of the order of other subleading corrections not included in our anal- 
ysis, and lead to a very small modification of the current Tevatron limits (a small shift, 
smaller than Atan/5 ~ 1, in the limit on tan/5). They have been neglected in the CDF 
analysis of the a{pp — > 0) x BR({/) r+r~) process. We shall omit these corrections in the 
analytical formulae presented in this section and also in the numerical analysis below. How- 
ever, including them would be straightforward (leading, as discussed above, to a very small 
modification of the allowed parameter space). Hence, the total production rate of bottom 
quarks and r pairs mediated by the production of a CP-odd Higgs boson in the large tan /? 
regime is approximately given by 

a(bbA) X BR(A hi) - a(bbA)sM ^ x . , (19) 

^ ^ ^ ^ ^ ^ il + A,f (l + A,)V9 ^ ^ 

aigg,bb^A)xBR{A^T+T-) ~ aigg,bb ^ A)sm , (20) 

(i + iXb) + y 

where a{bbA)sM and cr^gg, bb A)sm denote the values of the corresponding SM Higgs 
boson production cross sections for a Higgs boson mass equal to M4. 

As a consequence, the bb production rate depends sensitively on A;, because of the factor 
1/(1 -|- Afe)^, while this leading dependence on Af, cancels out in the r+r~ production rate. 
There is still a subdominant parametric dependence in the r+r" production rate on A;, that 
may lead to variations of a few tens of percent of the r-pair production rate (compared to 
variations of the rate by up to factors of a few in the case of bottom-quark pair production). 

The formulae above apply, within a good approximation, also to the non-standard CP- 
even Higgs boson in the large tan f3 regime. Indeed, unless the CP-odd Higgs mass is within 
a small regime of masses of about m™'^^ ~ 120-130 GeV, the mixing of the two CP-even 
Higgs bosons is small and, for Ma > m™^^ {Ma < m™'^^), cos a ^ sin/5 (sin a ~ — sin/5). 
In addition, this non-standard Higgs boson becomes degenerate in mass with the CP-odd 
Higgs scalar. Therefore, the production and decay rates of H (h) are governed by similar 
formulae as the ones presented above, leading to an approximate enhancement of a factor 
2 of the production rates with respect to the ones that would be obtained in the case of 
the single production of the CP-odd Higgs boson as given in eqs. (fT^. (|^. The same is 
true in the region where all three neutral Higgs bosons are approximately mass-degenerate, 
ruh — ruH — Ma- In this case the combined contribution of h and H to the production and 
decay rates approximately equals the contribution of the CP-odd Higgs boson. 

Besides the effects discussed above, additional radiative corrections can be important 
in the search for non-standard MSSM Higgs bosons. In particular, there are radiative cor- 
rections to the mass difference between mn (the non-SM like CP-even Higgs boson) and 
Ma [52]. If the two states are roughly mass-degenerate, one obtains a factor of 2 in the 
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production rate, as outlined above. In the case of small mixing between the two CP-even 
states, the mass difference is approximately given by 



m 



H 



Mi 



G p 



^SUSY 



+ 



ml 



(/^ At? 

^SUSY 



(21) 



where Gp is the Fermi constant, and denotes the running quark masses. For large values 
of ifiAt) and/or large values of (fiAb) and tan/?, the mass difference becomes so large that 
the signals arising from the production of the CP-even and the CP-odd Higgs bosons can 
no longer simply be added, leading to a modification of the Higgs search sensitivity at the 
Tevatron and the LHC A further important set of corrections are contributions from the 
sbottom sector giving rise to a large downward shift in the mass of the light CP-even Higgs 
boson, 

ft / P 1 (22) 



Sm. ~ 



AV27cHl + Aby VMsusY. 

For large values of n and tan/3 these corrections can shift the prediction for m/j below 
the experimental bound from LEP [1,2]. This can happen in particular for small mixing 
in the stop sector, for which the LEP bounds exclude a significant part of the parameter 
space. Finally, there are radiative corrections affecting the mixing of the two CP-even 
Higgs states that are not included in the above expressions. In particular, bottom- Yukawa 
induced corrections lead to an enhancement (suppression) of the mixing between the CP- 
even Higgs bosons for large and negative (positive) values of m^/(l + A;,)"^ fi^Ab/M^^^Y- 
enhancement of the mixing between h and H implies that the mass difference between mn 
and Ma is pushed to larger values than those given in eq. ()2H) and that mt receives a further 
downward shift in addition to the correction in eq. 



We now turn to the production and decay processes of the charged Higgs boson. In 
the MSSM, the masses and couplings of the charged Higgs boson in the large tan/? regime 
are closely related to the ones of the CP-odd Higgs boson. The tree-level relation M|^± = 
M\ + M^ receives sizable corrections for large values of tan /?, /x. At and Ab, 



Ml 



+ 



M^ 



3 Gp 



8^2 



2 2 

m^ mi 
(I + A5) 



4 log 



^ ) + 2Atb 



+ 



ml 



;l+A.)^ 



(23) 



/i 



SUSY 



with [18] 



A 



tb 



6/i2 



ifi'-AbAtf , 3{At + Ab 



M"^ 



M"^ 



+ 



^^-'SUSY 



(24) 



The coupling of the charged Higgs boson to a top and a bottom quark at large values of 
tan /3 is governed by the bottom Yukawa coupling and is therefore affected by the same A;, 
corrections that appear in the couplings of the non-standard neutral MSSM Higgs bosons [29]. 

The relevant channels for charged Higgs boson searches depend on its mass. For values 
of Mh± smaller than the top-quark mass, searches at hadron colliders concentrate on the 
possible emission of the charged Higgs boson from top-quark decays. In this case, for large 
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values of tan /?, the charged Higgs decays predominantly into a r lepton and a neutrino, i.e. 
one has to a good approximation 



BR(i/± ^ TVr) ~ 1 ■ (25) 

The partial decay width of the top quark into a charged Higgs and a bottom quark is propor- 
tional to the square of the bottom Yukawa coupling and therefore scales with tan^ /^/(l-l-Afe)^, 
see e.g. Ref. [29]. 

For values of the charged Higgs mass larger than mt, instead, the most efficient production 
channel is the one of a charged Higgs associated with a top quark (mediated, for instance, by 
gluon-bottom fusion). In this case, the production cross section is proportional to the square 
of the bottom-quark Yukawa coupling. The branching ratio of the charged Higgs decay into 
a T lepton and a neutrino is, apart from threshold corrections, governed by a similar formula 
as the branching ratio of the decay of the CT'-odd Higgs boson into r-pairs, namely 



where the factor (1 — r^)^ is associated with threshold corrections, and = m|/M^±. 



(26) 



As mentioned above, our numerical analysis will be based on the complete expressions for 
the Higgs couplings rather than on the simple approximation formulae given in this section. 



3 Interpretation of cross section limits in MSSM 
scenarios 

3.1 Limits at the Tevatron 

The DO and CDF Collaborations have recently published cross section hmits from the Higgs 

search at the Tevatron in the channel where at least three bottom quarks are identified 
in the final state (660, bb) [7] and in the inclusive channel with r+r~ final states 
{pp — > — > r+T") [8]. The CDF Collaboration has also done analyses searching for a 
charged Higgs boson in top-quark decays [9]. For a Higgs boson with a mass of about 
120 GeV, the DO Collaboration excludes a cross section of about 30 pb in the 660, — > 66 
channel with a luminosity of 260 pb~^ [7], and the CDF Collaboration excludes a cross 
section of about 15 pb in the pp ^ (p ^ r^r^ channel with a luminosity of 310 pb^^ [8]. 
While the cross section for a SM Higgs boson is significantly below the above limits, a large 
enhancement of these cross sections is possible in the MSSM. 

It is therefore of interest to interpret the cross section limits within the MSSM parameter 
space. Since the Higgs sector of the MSSM is characterised by two new parameters at lowest 
order, conventionally chosen as Ma and tan/9, one usually displays the limits in the M^- 
tan/3 plane (for CP- violating scenarios one normally chooses the M^-tan/? plane). As the 
whole structure of the MSSM enters via radiative corrections, the limits in the M^-tan/? 
plane depend on the other parameters of the model. One usually chooses certain benchmark 
scenarios to fix the other MSSM parameters [3,4]. In order to understand the physical 



9 



meaning of the exclusion bounds in the M^-tan (3 plane it is important to investigate how 
sensitively they depend on the values of the other MSSM parameters, i.e. on the choice of 
the benchmark scenarios. 



3.1.1 Limits from the process bb 

The DO Collaboration has presented the limits in the M^-tan f3 plane obtained from the 
bb(j), (f) bb channel for the m™*^^ and no-mixing scenarios as defined in Ref. [3]. The m™*^^ 
scenario according to the definition of Ref. [3] reads 



mt = 174.3 GeV, 
MsusY = 1000 GeV, 
/i = -200 GeV, 
Ms = 200 GeV, 

= 2 MsusY (FD calculation), 
Xf^ = VQMsvsy (RG calculation) 
Ab = At, 

m-g = O.SMsusY . (27) 

The no-mixing scenario defined in Ref. [3] differs from the m'^^^ scenario only in 

Xt = (FD/RG calculation) . (28) 

The condition Af, = At implies that the different mixing in the stop sector gives rise to a 
difference between the two scenarios also in the sbottom sector. The definition of the m™'^^ 
and no-mixing scenarios given in Ref. [3] was later updated in Ref. [4], see the discussion 
below. 

For their analysis, the DO Collaboration has used the following approximate formula [7], 

a(bb(P) X BR(0 ^bb) = 2 a(bb^)sM ^ x . , (29) 

^ ' ' ^ ^^'"^ (1 + Ab)' {l + Abf + 9 ^ ' 

which follows from eq. (fT^ and the discussion in Sect. 12.31 The cross section (j{bb(t>)sM has 
been evaluated with the code of Ref. [48], while A;, has been calculated using CPsuperH [38]. 
From the discussion in Sect. 12.31 it follows that the choice of negative values of fi leads to an 
enhancement of the bottom Yukawa coupling and therefore to an enhancement of the signal 
cross section in eq. ()29p. For tan (3 = 50 the quantity A;, takes on the following values in the 
m™^^ and no-mixing scenarios as defined in eqs. ()27|) . ()28|) . 

mj^^ scenario, = -200 GeV, tan /? = 50 : A^ = -0.21 , (30) 
no-mixing seen., = -200 GeV, MsusY = 1000 GeV, tan /3 = 50 : A^ = -0.10 . (31) 

While the 0{as) contribution to A;,, see eq. fll5|) . is practically the same in the two scenarios, 
the 0{at) contribution to Af, in the m™^^ scenario differs significantly from the one in the no- 
mixing scenario. In the m™'^^ scenario the 0{at) contribution to A;, is about as large as the 
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0{as) contribution. In the no-mixing scenario, on the other hand, the 0{at) contribution to 
Af, is very small, because At is close to zero in this case. Reversing the sign of /i in eqs. (jHUj) . 
(jnH) reverses the sign of A;,, leading therefore to a significant suppression of the signal cross 
section in eq. (j29|l for the same values of the other MSSM parameters. 

The predictions for bb(f), (p bb evaluated with FeynHiggs have been compared with the 
exclusion bound for a x BR as given in Ref. [7]. As mentioned above, in our analysis we 
use the full Higgs couplings obtained with FeynHiggs rather than the approximate formula 
given in eq. ()29p. Similar results would be obtained with CPsuperH. 

The impact on the limits in the M^-tan (3 plane from varying fi while keeping all other 
parameters fixed can easily be read off from eq. For a given value of the CP-odd mass 

and tan/3, the bound on a{bb(f)) x BR(0 bb) provides an upper bound on the bottom- 
quark Yukawa coupling. The main effect therefore is that as n varies, the bound on tan (3 
also changes in such a way that the value of the bottom Yukawa coupling at the boundary 
line in the M^-tan (3 plane remains the same. 

The dependence of the limits in the M^-tan f3 plane obtained from the process bb(f), (p bb 
on the parameter fi is shown in Fig. ^ The limits for /i = —200 GeV in the m™^^ and no- 
mixing scenarios, corresponding to the limits presented by the DO Collaboration in Ref. [7], 
are compared with the limits arising for different /i values, n = -1-200, ±500, ±1000 GeV. 
Fig. ^illustrates that the effect of changing the sign of fi on the limits in the M^-tan f3 plane 
obtained from the process bb(j), (p —>■ bb is quite dramatic. In the m^^^ scenario the exclusion 
bound degrades from about tan (3 = 50 for Ma = 90 GeV in the case of /i = —200 GeV to 
about tan /5 = 90 for Ma = 90 GeV in the case of /i = ±200 GeV. We extend our plots to 
values of tan j3 much larger than 50 mainly for illustration purposes; the region tan (3 ^ 50 
in the MSSM is theoretically disfavoured, if one demands that the values of the bottom 
and r Yukawa couplings remain in the perturbative regime up to energies of the order of 
the unification scale. The situation for the bottom- Yukawa coupling can be ameliorated for 
large positive values of /i due to the A;, corrections. The curves for /i = ±500, ±1000 GeV 
do not appear in the plot for the m™*^^ scenario, since for these /i values there is no tan/3 
exclusion below tan/3 = 130 for any value of Ma- On the other hand, the large negative 
values of fi shown in Fig. ^ fi = —500, —1000 GeV, lead to an even stronger enhancement 
of the signal cross section than for fi = —200 GeV and, accordingly, to an improved reach 
in tan/3. It should be noted that for fi = —500, —1000 GeV the bottom Yukawa coupling 
becomes so large for tan/3 ^ 50 that a perturbative treatment would no longer be reliable 
in this region. 

In the no-mixing scenario, where the absolute value of A;, is smaller, the exclusion bound 
is shifted from about tan /3 = 55 for Ma = 90 GeV in the case of = —200 GeV to about 
tan p = 75 for Ma = 90 GeV in the case of = ±200 GeV. For /x = ±500 GeV, no 
excluded region can be estabhshed below tan/3 = 100. As above, large negative values of 
/i, i.e. /i = —500, —1000 GeV, result in an improved reach in tan/3 compared to the value 
/i = —200 GeV chosen by the DO Collaboration. 

The variation with the sign and absolute value of the /x-parameter shows the strong 
dependence of the limits in the M^-tan (3 plane on the strength of the bottom-quark Yukawa 
coupling and hence on the supersymmetric parameter space. The sensitive dependence of the 
process bbcp, cp —> bb on the bottom-quark Yukawa coupling is not specific to the particular 
benchmark scenarios considered here. Keeping the ratio of /im^/MgugY fixed but varying /i 
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Figure 1: Change in the hmits obtained from the bb(j), (p ^ bb channel in the m™^^ (left) and 
no-mixing (right) benchmark scenarios for different values of /i. The value fi = —200 GeV 
was chosen by the DO Collaboration in Ref. [7]. The other curves indicate the corre- 
sponding limfts for fi = +200, ±500, ±1000 GeV. The curves for /i = ±500, ±1000 GeV 
(yU = ±1000 GeV) do not appear in the left (right) plot for the m™^^ (no-mixing) scenario, 
since for these /x values there is no tan (3 exclusion below tan (3 = 130 for any value of Ma- 

and rrig independently will lead to similar results as those shown here. A scenario where large 
compensations are possible between the two contributions entering A^, see eq. flT^ . will be 
discussed below. Scenarios with different values of the other supersymmetric parameters 
(besides the ones entering A;,) will reproduce a similar behaviour as those discussed here. 

In Ref. [4] the definition of the m™'^^ and no-mixing scenarios given in Ref. [3] has been 
updated, and the "small ctes" scenario and the "gluophobic Higgs scenario" have been pro- 
posed as additional scenarios for the search for the light CP-even Higgs boson at the Teva- 
tron and the LHC. The sign of fi in the m™'^^ and no-mixing scenarios has been reversed to 
fi = ±200 GeV in Ref. [4]. This leads typically to a better agreement with the constraints 
from {g — 2)^. Furthermore, the value of Msusy in the no-mixing scenario was increased 
from 1000 GeV [3] to 2000 GeV in order to ensure that most of the parameter space of this 
scenario is in accordance with the LEP exclusion bounds [1,2]. 

Another scenario defined in Ref. [4] is the "constrained-m™'^^" scenario. It differs from 
the m™^'' scenario as specified in Ref. [4] by the reversed sign of Xt, 

—2 Msusy (FD calculation), 
—VQ MsvsY (RG calculation), 

±200 GeV . (32) 

For small Ma and minimal flavor violation this results in better agreement with the con- 
straints from BR(6 — > 57). For large tan (3 one has At ^ Xt, thus At and rrig have opposite 
signs. This can lead to cancellations in the two contributions entering A;,, see eq. f|T^ . In 
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contrast to the mf^^ scenario, where the two contributions entering A;, add up, see eq. (IHUj). 
the constrained-m™^^ scenario typically yields relatively small values of A;, and therefore 
a correspondingly smaller effect on the relation between the bottom-quark mass and the 
bottom Yukawa coupling, e.g. 



constrained-mjf^'' scenario, /i = +200 GeV, tan /? = 50 : A^ 



-0.001 



(33) 



For large values of the compensations between the two terms entering A;, are less efficient, 
since the function / in the second term of eq. ()3()|1 scales like l//i^ for large 

We now study the impact of the benchmark definitions of Ref. [4] on the limits in the 
M^-tan (3 plane arising from the 660, (p bb channel. The left plot in Fig. El shows the 
effect of changing Msusy = 1000 GeV to Msusy = 2000 GeV in the no-mixing scenario for 
fi = ±200 GeV. Due to the heavier scalar bottoms in the case of Msusy = 2000 GeV the 
effect of the Af, corrections is suppressed as compared to the benchmark definition in Ref. [3]. 
This leads to a shift of the limits in the M^-tan (3 plane of about A tan (3 = 5-10 for a given 
value of Ma- The right plot of Fig. |21 shows for Msusy = 2000 GeV the variation of the limits 
with yu. In this case even for /i = +1000 GeV a tan (3 exclusion limit can be established below 
tan/3 = 130, in contrast to the scenario with Msusy = 1000 GeV, see Fig. [H 
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Figure 2: Variation of the limits obtained from the bb(f), (j) bb channel in the no- 
mixing scenario for different values of Msusy and /x. The left plot shows the results for 
Msusy = 1000, 2000 GeV and = ±200 GeV, while in the right plot the results for 
Msusy = 2000 GeV and /i = ±200, ±500, ±1000 GeV are given. 

The results in the constrained-m™'^^ scenario are displayed in Fig. El (left). The results 
are shown for n = ±200, ±500 GeV. As expected from the discussion above, the obtained 
limits are relatively stable against the variation of fi. For fi = +500(— 500) GeV the tan/? 
limit is significantly weaker (stronger) than for smaller values of as a consequence of the 
less efficient cancellation of the two contributions to A^ discussed above. Nevertheless, the 
limits obtained for |//| < 500 GeV are weaker than those for the m™'^'' scenario with negative 
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H, but stronger than those for positive /i. The curve for n = +1000 GeV is not shown in 
the plot, since for this value there is no tan/5 exclusion below tan/5 = 130 for any value of 
Ma. For fi = —1000 GeV, on the other hand, the radiative corrections lead to a large mass 
splitting between the CP-odd and CP-even Higgs boson masses so that the approximation 
of adding the two signal cross sections is no longer valid, see the discussion in Sect. 12.31 A 
more detailed study would be necessary to incorporate also the case of larger Higgs boson 
mass splittings. 



130 

120 

110 

100 

90 

80 

CO. 70 
c 
to 



60i 
50 
40 
30 
20 
10 

90 



constr. , Tevatron: bbc^, ^ -> bb 
. n = -500 GeV 



: -200 GeV 
: -1-200 GeV 
: -1-500 GeV 



CO. 

c 
03 



130 
120 
110 
100 
90 
80 
70 
60 
50 
40 
30 
20 
10 



100 



110 120 

[GeV] 



130 



140 



150 



_l 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 








constr. m^™*^, Tevatron: ()> -> tt, full A^^ : 






. 11 = -500 GeV 






- = -200 GeV J 






- n = -1-200 GeV J 






H = -1-500 GeV : 






- n = -1-1 000 GeV ^ 


i, , , , 1 , , 




, 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ,, 



90 100 110 120 130 140 150 160 

[GeV] 



170 180 190 200 



Figure 3: Left: Variation of the limits obtained from the bb(j), (f) ^ bb channel in the 
constrained-m™^^ scenario for different values of fi. Right: Variation of the limits obtained 
from the pp ^ (f) t^t" channel in the constrained-m™'^^ scenario for different values of /i. 



3.1.2 Limits from the process pp ^ (j) t^t 

The limits obtained from the pp ^ (f) ^ r^r^ channel by the CDF Collaboration were 
presented in the M4-tan (3 plane for the m™'^^ and no-mixing scenarios as defined in Ref. [4] 
and employing two values of the /i parameter, fi = ±200 GeV. According to the discussion 
in Sect. 12. 3| the limits obtained from the pp —* (p ^ t^t^ channel are expected to show 
a weaker dependence on the sign and absolute value of than the limits arising from the 
660, — > 66 channel. On the other hand, for large values of tan j3 and negative values of 
/i, the large corrections to the bottom Yukawa coupling discussed above can invalidate a 
perturbative treatment for this channel. 

The MSSM prediction for a{pp — 0) x BR(0 t'^t~) as a function of tan/5 has 
been evaluated by the CDF collaboration using the HIGLU program [55] for the gluon fusion 
channel. The prediction for the 66 — >■ + X channel was obtained from the NNLO result in 
the SM from Ref. [50], and [a x BR]j^ggp^ / W >^ ^^]sm calculated with the FeynHiggs 
program [34-37]. While the full correction to the bottom Yukawa correction was taken 
into account in the 66 ^ + X production channel and the t^t~ branching ratios. 
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the public version of the HIGLU program [55] does not include the A;, correction for the 
bottom Yukawa coupling entering the bottom loop contribution to the gg ^ (p production 
process. In order to treat the two contributing production processes in a uniform way, the 
Ab correction should be included (taking into account the 0{as) and the 0{at) parts, see 
eq. (fTH|) ) in the gg (p production process calculation. For the large value of Msusy chosen 
in the m™'^^ and no-mixing benchmark scenarios other higher-order contributions involving 
sbottoms and stops can be neglected (these effects are small provided Msusy ^ 500 GeV). 

We therefore begin the investigation of the pp ^ (p t^t^ channel by analyzing the 
impact of including or omitting the A;, correction in the gg ^ cp production process. In 
order to get a qualitative understanding of the variation of the limits on tan (3 for a given 
value of Ma induced by the inclusion of the A;, corrections in the gluon fusion channel, it 
is again useful to employ the simple approximate formulae given in Sect. 12.31 As discussed 
above, the production cross sections may be approximately obtained from the SM ones by 
including a simple rescaling by a factor tan^/3/(l + A;,)^. Hence, defining ab and ag as the 
SM cross sections for the 6-quark associated and gluon fusion production of Higgs bosons, 
respectively, we get 

a(pp-0)xBR(0^r+r-)|^^^^^^ - + ^ (TT^^ ' ^'^^ 
a(pp-.0)xBR(0^r+r-) ^ U ^ {I + Abf) x , (35) 



partial A, V " ^ ^ (1 + Ab)^ + 9 

where "full A;," denotes the case where the A;, correction is incorporated in both the h- 
quark associated and the gluon fusion production processes (and the (p — r^r~ branching 
ratio) , while "partial Ab" denotes the case where the A^ correction is omitted in the gg ^ cp 
production process. The expression above shows that, for positive values of /x, for which 
Ab > 0, the omission of the correction to the gg ^ cp process leads to an enhancement of 
the total production cross section with respect to the value obtained when these corrections 
are included. For negative values of /i, instead, the situation is reversed. 

The cross section for the gg ^ cp production process including the Ab correction can be 
obtained by a simple rescaling of the HIGLU result for the SM. The cross section for the SM 
production rate involving the 6-quark loop alone is rescaled with V{h bb)MSSM/^{h — > 
bb)sM, where Ab enters the calculation of T{h — >• 66)mssm (SM-QCD corrections to the b- 
quark mass factorize and drop out in this ratio). As stressed above, it has been checked that 
the t-quark loop gives only a negligible contribution in the MSSM for tan f3 ^ 20. The loops 
involving scalar tops and bottoms, beyond those included in the Ab corrections, also give 
small contributions due the heavy scalar masses in the benchmark scenarios. 

The comparison of the "partial Ab" and the "full Ab" results is shown in Fig. |3 The 
evaluation of a x BR in the MSSM has been performed by using FeynHiggs for rescaling the 
HIGLU results for the gg ^ cp production cross sections and the SM results for the bb cp + X 
channel [50] by the appropriate MSSM correction factor, as outlined above. The impact of 
the additional contribution can be read off from Fig. 0] by comparing the results where the 
Ab corrections are omitted in the gg (p production cross sections ( "partial Ab" ) with the 
results where the Ab corrections have been taken into account everywhere in the production 
and decay processes ("full Ab"). The effect on the exclusion bounds in the M^-tan/3 plane is 
seen to be quite significant. While in the case where the Ab corrections are neglected in the 
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gg ^ (j) production cross sections the strongest exclusion bounds are obtained for positive 
values of /x, the inclusion of the A;, corrections reverses this situation. As explained above, 
the inclusion of the A;, corrections to the gg ^ (p production process leads to a larger cross 
section (and correspondingly to a stronger tan P bound) in the case of negative /i, while the 
cross section is suppressed for positive values of /i. The corresponding shifts of the exclusion 
limits in the M^-tan /3 plane amount up to A tan /5 ~ 10 for the m™*^^ scenario. In the no- 
mixing sceanrio (defined according to Ref. [4]) the effect is less pronounced because of the 
smaller numerical value of A;,, giving rise to shifts in the exclusion limits up to Atan/3 ~ 5. 




[GeV] [GeV] 



Figure 4: Impact of including or omitting the A^ correction in the gg ^ cj) production 
process on the limits obtained from the pp ^ (p, (p ^ t^t~ channel. The results are shown 
for /i = ±200 GeV in the m™^^ (left) and no-mixing (right) benchmark scenarios [4]. 

Following our analysis, the CDF Collaboration has adopted the prescription outlined 
above for incorporating the A;, correction into the gg ^ (p production process. The limits 
given in Ref. [8] are based on the MSSM prediction where the A;, correction is included 
everywhere in the production and decay processes (see e.g. Ref. [56] for a previous analysis). 

We next turn to the discussion of the sensitivity of the limits obtained from the pp 
(p T^T^ channel (including the Af, correction in all production and decay processes) on 
the sign and absolute value of ^. As discussed above, similar variations in the exclusion 
limits will occur if the absolute values of /i, m^. At and Msusy are varied, while keeping the 
ratios appearing in A;, constant. The results are given in Fig. |Slfor the rnf^'^^ scenario (left) 
and the no-mixing scenario (right). In the rn^^^ scenario we find a sizable dependence of the 
tan [3 bounds on the sign and absolute value of ^i? The effect grows with Ma and, for the 
range of parameters explored in Fig. leads to a variation of the tan (3 bound larger than 

^For /i = —300 GeV the curve stops at around tan/3 = 75 because the bottom Yukawa coupUng becomes 
very large, leading to instabilities in the calculation of the Higgs properties. For the same reason, even more 
negative values of n are not considered here. 
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A tan j3 ~ 30. In the no-mixing scenario the effect is again smaller, but it can still lead to a 
variation of the tan (3 bounds by as much as Atan/5 ~ 10. 




[GeV] [GeV] 



Figure 5: Variation of the limits obtained from the pp ^ (f) ^ r^r channel at the Tevatron 
in the m™^^ (left) and no-mixing (right) benchmark scenarios for different values of /i. 

The results obtained in the constrained-m™^^ scenario are shown in Fig. El (right). As 
expected, the exclusion limits in this scenario are very robust with respect to varying fi. All 
values of /i result practically in the same tan (3 exclusion bounds. The lines not visible in the 
plot are actually covered by a line of another /j, value. For /i = —1000 GeV, the radiative 
corrections lead to a large mass splitting between the CP-odd and CP-even Higgs boson 
masses, see the discussion above. 

3.1.3 Limits from the process pp ^ ti ^ H'^W^ bb, rur 

For the charged Higgs search channel at the Tevatron [9], pp tt ^ H'^W^bb, — * rur 
the variation of the cross section with fi is driven by the impact of the Af, correction on 
BR(t HH) [29]. The decay width r(t H^b) is proportional to tan^ /3/(l + Afe)^, leading 
to an expression for the branching ratio in analogy to eq. (jl7|) . Accordingly, a positive A;, 
leads to a suppression of BR(t H^b), while a negative A;, leads to an enhancement. 

For a fixed value of Mh±, the value of Ma is driven to rather small values because of 
the tree-level relation M\ = Mfj± — M^. For large values of tan/3. At and Af, this effect 
is further enhanced by the higher-order corrections in eq. (|23|). Consequently, in the region 
of small Mh± and large tan/3 currently probed at the Tevatron [9] the corresponding M4 
values tend to be further reduced with respect to the already small tree-level values^ and 
are in general below the LEP exclusion bound [2]. Therefore, this channel at present is less 

•^This effect is avoided if the parameters At, and ^ are such that Atb in eq. (|24|l becomes sufficiently 
large and negative. This can be realized, for instance, for values of Mh± > 120 GeV and At = —Ab in the 
constrained-m™'*'^ scenario, for small values of fi. 
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relevant for obtaining exclusion limits in the M^-tan (3 plane than the neutral Higgs-boson 
search channels discussed above. It is expected to become more competitive, however, with 
increasing luminosity collected in Run II of the Tevatron. 

3.2 Prospects for Higgs sensitivities at the LHC 

The most sensitive channels for detecting heavy MSSM Higgs bosons at the LHC are the 
channel pp H/A + X, H/A t^t~ (making use of different decay modes of the two r 
leptons) and the channel tH^,H'^ rUr (for Mi{± > rrit) [57,58]. We consider here the 
parameter region Ma 3> M^, for which the heavy states H , A are widely separated in mass 
from the light CP-even Higgs boson h. Here and in the following we do not discuss search 
channels where the heavy Higgs bosons decay into supersymmetric particles, which depend 
very sensitively on the model parameters [58-60], but we will comment below on how these 
decays can affect the searches with bottom-quarks and r-leptons in the final state. 

3.2.1 Discovery region for the process pp — > H/A + X, H/A — > r+r" 

To be specific, we concentrate in this section on the analysis carried out by the CMS Col- 
laboration [58,61]. Similar results for this channel have also been obtained by the ATLAS 
Collaboration [57, 62] . In order to rescale the SM cross sections and branching ratios, the 
CMS Collaboration has used for the branching ratios the HDECAY program [63] and for the 
production cross sections the HIGLU program [55] {gg H/A) and the HQQ program [64] 
{gg bbH). In the HDECAY program the corrections are partially included for the decays 
of the neutral Higgs bosons (only the 0{as) contribution to is included, see eq. (HH)))- 
The HIGLU program (see also the discussion in Sect. I3.1.2|) and HQQ, on the other hand, do 
not take into account the corrections to the bottom Yukawa coupling.^ The prospective 5a 
discovery contours for CMS (corresponding to the upper bound of the LHC "wedge" re- 
gion, where only the light CP-even Higgs boson may be observed at the LHC) have been 
presented in Refs. [58,61] in the M^-tan/^ plane, for an integrated luminosity of 30 fb~^ 
and 60 fb~^. The results were presented in the m™*^^ scenario and for different fi val- 
ues, = —200, -1-300, -1-500 GeV. It should be noted that decays of heavy Higgs bosons 
into charginos and neutralinos open up for small enough values of the soft supersymmetry- 
breaking parameters M2 and /x. Indeed, the results presented in Refs. [58,61] show a degra- 
dation of the discovery reach in the M^-tan (3 plane for smaller absolute values of /z, which is 
due to an enhanced branching ratio of H, A into supersymmetric particles, and accordingly 
a reduced branching ratio into r pairs. 

We shall now study the impact of including the Ab corrections into the production cross 
sections and branching ratios for different values of fi. The inclusion of the A;, corrections 
leads to a modification of the dependence of the production cross section on tan/3, as well 
as of the branching ratios of the Higgs boson decays into r+r^. For a fixed value of M^, 
the results obtained by the CMS Collaboration for the discovery region in tan/5 can be 
interpreted in terms of a cross section limit using the approximation of rescaling the SM rate 

''Since HQQ is a leading-order program, non-negligible changes can also be expected from SM-QCD type 
higher-order corrections. 
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for the pp ^ H + X, H ^ t^t process by the factor 

^ 2. BR(iJ^r+r-)cMS + BR(A^r+r-)cMS 

tan /3cMS X — — -— . (36) 

BR(ii T+T )sM 

In the above, tan /?cms refers to the value of tan (3 on the discovery contour (for a given 
value of Ma) that was obtained in the analysis of the CMS Collaboration with 30 fb^^ [58]. 
These tan [3 values as a function of Ma correspond to the edge of the area in the M^-tan [3 
plane in which the signal pp — > H/A + X, H/A t^t~ is visible (i.e. the upper bound of 
the LHC wedge region). The branching ratios BR{H t~^t~)cms and BR(74 — > t~^t~)cms 
in the CMS analysis have been evaluated with HDECAY, incorporating therefore only the 
gluino-sbottom contribution to A;,. 

After including all A;, corrections, we evaluate the pp H/A + X, H/A t^t~ process 
by rescaling the SM rate with the new factor, 

tan^/3 BR(ij' ^ r+r') + BR(A ^ r+r') 
(1 + Afe)2 BR{H ^ r+r-)sM ' ^ ' 

where A^ depends on tan/3. The quantities have been evaluted with FeynHiggs, allowing 
also decays into supersymmetric particles. The resulting shift in the discovery reach for the 
pp — > H/A + X, H/A T^T^ channel can be obtained by demanding that eq. (jHUj) and 
eq. pTj) should give the same numerical result for a given value of Ma- 

This procedure has been carried out in two benchmark scenarios for various values of /x. 
The results are shown in Fig. IHlfor the m™'^^ scenario (left) and for the no-mixing scenario 
(right). The comparison of these results with the ones obtained by the CMS Collabora- 
tion [58,61] shows that for positive values of ^ the inclusion of the supersymmetric radiative 
corrections leads to a slight shift of the discovery region towards higher values of tan i.e. 
to a small increase of the LHC wedge region. For = —200 GeV the result remains ap- 
proximately the same as the one obtained by the CMS Collaboration. Due to the smaller 
considered tan/? values compared to the analysis of the Tevatron limits in Sect. 13. iT^ the 
corrections to the bottom Yukawa coupling from Af, are smaller, leading to a better pertur- 
bative behavior. As a consequence, also the curves for /i = —500, —1000 GeV are shown in 
Fig. El 

The change in the upper limit of the LHC wedge region due to the variation of /i does 
not exceed A tan/? ~ 8. As explained above, this is a consequence of cancellations of the 
leading Af, effects in the Higgs production and the Higgs decay. Besides the residual A^ 
corrections, a further variation of the bounds is caused by the decays of the heavy Higgs 
bosons into supersymmetric particles. For a given value of /x, the rates of these decay modes 
are strongly dependent on the particular values of the weak gaugino mass parameters M2 
and Ml. In our analysis, we have taken M2 = 200 GeV, as established by the benchmark 
scenarios defined in Ref. [4], while Mi ~ 100 GeV. Since the Higgs couplings to neutralinos 
and charginos depend strongly on the admixture between higgsino and gaugino states, the 
rate of these processes is strongly suppressed for large values of > 500 GeV. In general, 
the effects of the decays H/A x'iX^^xtxT o^ily play a role for M4 > + Mi. Outside 
this range the cancellations of the A;, effects result in a very weak dependence of the rates 
on /i. 



19 




200 



300 



I I I I I I I I I I 
400 500 

[GeV] 



600 



700 



800 



200 



300 



I I I I I I I I I I 
400 500 

[GeV] 



600 



700 



800 



Figure 6: Variation of the 5a discovery potential for the pp H/A + X, H/A t^t~ 
process at the LHC in the mf^^^ (left) and no-mixing (right) benchmark scenarios for different 
values of /i. 

The combination of the effects from supersymmetric radiative corrections and decay 
modes into supersymmetric particles gives rise to a rather complicated dependence of the 
discovery contour on /i. This is illustrated in Fig. [7| (left), where the discovery contour for 
the pp H/A + X, H/A —>■ t^t~ process is shown as a function of /i in the m™^^ scenario 
for different values of Ma- As explained above, for Ma ^ |/^| + the decay modes into 
supersymmetric particles have a significant impact, while outside this region the dependence 
on IX is rather weak. 

In Fig. ISl (left) we show the results for the constrained-m™'^^ scenario, see eq. ()32|) . The 
variation of the discovery contour in the M^-tan/5 plane with /i is completely driven in 
this case by the additional decay channels of the heavy Higgs bosons into charginos and 
neutralinos. Correspondingly, the weakest sensitivity is obtained for the smallest values of 
^i = ±200 GeV. 

3.2.2 Discovery region for the process tH^^H^ rur 

For this process we also refer to the analysis carried out by the CMS Collaboration [58,65]. 
The corresponding analyses of the ATLAS Collaboration can be found in Refs. [57,66]. The 
results of the CMS Collaboration were given for an integrated luminosity of 30 fb~^ in the 
M4-tan/3 plane using the m™^^ scenario with ^ = —200 GeV. No A;, corrections were 
included in the gb —>■ tH^ production process [67] and the H^ — >• rvr decay [63]. 

In Fig. El we investigate the impact of including the Af, corrections into the production 
and decay processes and of varying /i. In order to rescale the original result for the discovery 
reach in tan/3 we have first evaluated the tan/3 dependence of the production and decay 
processes. If no supersymmetric radiative corrections are included, for a fixed Ma value, the 
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Figure 7: Variation of the 5a discovery contours as a function of the parameter // in the m™^ 
scenario for the pp — > H/A + X, H/A — > r+r" process (left) and the H^ — > tu^ process 
(right). 




Figure 8: Variation of the 5a discovery contours for different values of // in the constrained- 
j^max scenario for the pp — > H/A + X, H/A — > t+t~ process (left) and the H^ — > ruj. 
process (right). 
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Figure 9: Variation of the 5a discovery contours obtained from the tH^, —>■ rv^- channel 
in the m™^^ (left) and no-mixing (right) benchmark scenarios for different values of ^. 

discovery potential can be inferred by using a rate approximately proportional to 

tan^ /5cMS X BR{H^ ^ tz/^cms • (38) 



Here tan/3cMS is given by the edge of the area in the M^-tan/^ plane in which the signal 
if^ — ^ TVr is visible, as obtained in the CMS analysis. The BR(if^ — ^ Ti'r)cMS has been 
evaluated with HDECAY. 

The rescaled result for the discovery contour, including all relevant A;, corrections, is 
obtained by demanding that the contribution 

X BR{H^ Tv;) , (39) 

where A;, depends on tan/?, is numerically equal to the one of eq. (IHHj). The quantities in 
eq. have been evaluated with FeynHiggs. 

This procedure has been carried out in two benchmark scenarios for various values of yU. 
The results are shown in Fig. IHlfor the m™^^ scenario (left) and for the no-mixing scenario 
(right). As a consequence of the cancellations of the leading A^ effects in the Higgs production 
and the Higgs decay the change in the discovery contour due to the variation of /i does not 
exceed A tan /5 ~ 10(6) in the m™'^^ (no- mixing) scenario. Also in this case there is a variation 
of the contour caused by decays into supersymmetric particles that, as in the neutral Higgs 
boson case, are only relevant for small values of |yu|. For completeness, in Fig. [7| (right) we 
show the corresponding variation of the discovery contour for the m™^^ scenario as a function 
of /i, for different values of Ma- Outside the region where the decays into supersymmetric 
particles are relevant the dependence on /x is relatively weak, but somewhat more pronounced 
than for the case of the neutral Higgs bosons H and A. 
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In Fig. IHl (right) we show the results for the constrained-m™^^ scenario, see eq. (j32j) . As in 
the case of the neutral Higgs bosons, the variation of the discovery contour in the M^-tan (3 
plane is completely driven by the additional decay channels of the heavy Higgs bosons into 
charginos and neutralinos. 



4 Benchmark Scenarios 



The benchmark scenarios defined in Ref. [4], which were mainly designed for the search for 
the light CP-even Higgs boson h in the CP-conserving case, are also useful in the search for 
the heavy MSSM Higgs bosons H, A and H^. In order to take into account the dependence 
on /i, which as explained above is particularly pronounced for the bb(j), cf) hh channel, we 
suggest to extend the definition of the m™^^ and no-mixing scenarios given in Ref. [4] by 
several discrete values of /i. The scenarios defined in Ref. [4] read 
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The constrained-m™^^ scenario differs from eq. (j4Uj) only by the reversed sign of Xf. While the 
positive sign of the product (/i M2) results in general in better agreement with the [g — 2)^ 
experimental results, the negative sign of the product (fJ^Af) yields in general (assuming 
minimal flavor violation) better agreement with the BR(6 57) measurements. 

Motivated by the analysis in Sect. El we suggest to investigate the following values of /i 

fx = ±200, ±500, ±1000 GeV , (43) 

allowing both an enhancement and a suppression of the bottom Yukawa coupling and taking 
into account the limits from direct searches for charginos at LEP [68]. As discussed above, the 
results in the constrained-m™^^ scenario are expected to yield more robust bounds against 
the variation of /i than in the other scenarios. It should be noted that the values fi = 
—500, —1000 GeV can lead to such a large enhancement of the bottom Yukawa coupling 
that a perturbative treatment is no longer possible in the region of very large values of tan j3. 
Some care is therefore necessary to assess up to which values of fi reliable results can be 
obtained, see e.g. the discussion of Fig. El 

The value of the top-quark mass in Ref. [4] was chosen according to the experimental 
central value at that time. We propose to substitute this value with the most up-to-date 
experimental central value for m^. 

5 Conclusions 

In this paper we have analyzed the impact of supersymmetric radiative corrections on the 
current MSSM Higgs boson exclusion limits at the Tevatron and the prospective discovery 
reach at the LHC In particular, we have studied the variation of the exclusion and discovery 
contours obtained in different MSSM benchmark scenarios under changes of the higgsino 
mass parameter /i and the supersymmetry breaking parameters associated with the third 
generation squarks. These parameters determine the most important supersymmetric radia- 
tive corrections in the large tan/? region that are associated with a change of the effective 
Yukawa couplings of the bottom quarks to the Higgs fields (since the squarks are relatively 
heavy in the considered benchmark scenarios, other squark-loop effects are sub- dominant). 
These corrections had been ignored or only partially considered in some of the previous 
analyses of Higgs searches at hadron colliders. We have shown that their inclusion leads to 
a significant modification of the discovery and exclusion regions. 

We have investigated the exclusion bounds obtained from the Tevatron searches for non 
SM-like Higgs bosons in different channels. For the bb(j), cf) hh channel (0 = h, H, A) we 
find that the effects of the supersymmetric radiative corrections on the exclusion bounds in 
the M^-tan f3 plane are quite dramatic. While in the m™'^^ scenario the current data allow to 
rule out values of tan /? > 50 (35) for ^ 100 GeV if the higgsino mass parameter is chosen 
as /i = —200 GeV (—1000 GeV), hardly any bound on tan/? can be set if positive values 
of fi are chosen. The shifts are smaller, but still important, for the no-mixing benchmark 
scenario. We have shown that the constrained-m™*^^ scenario yields results that are much 
more stable against variations of fi than the other benchmark scenarios. 

For the inclusive channel with r"*"r~ final states, pp ^ (p ^ t~^t~ , compensations occur 
between large corrections to Higgs production and decay, so that the limits in the M^-tan /? 
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plane obtained from this channel turn out to be less affected by varying fj, than the ones 
from the associated production with bottom quarks. Nevertheless we have found that the 
exclusion limit is shifted by up to A tan /3 = 30 as a consequence of choosing different input 
values for We have investigated the impact of including the dominant supersymmetric 
radiative corrections to the gluon fusion production process, which had previously been 
omitted. The inclusion of these corrections leads to a shift of up to A tan/? = 10 in the 
exclusion limit. Following our analysis, the CDF Collaboration has adopted the prescription 
outlined in this paper for incorporating the correction into the gg ^ (p production process. 
The Tevatron experiments are expected to collect further data at higher luminosities, up to 
4-8 fb~^, in the next few years. This will extend the Tevatron MSSM Higgs boson discovery 
and exclusion reach in the Af^-tan (3 plane to lower values of tan decreasing the sensitivity 
of the obtained bounds to variations of the low energy supersymmetry mass parameters. 

For the LHC we have analyzed the channels pp — > H/A + X, H/A — > r+r~ and 
tH^, — > TUr, which yield the best sensitivities in the search for heavy MSSM Higgs 
bosons. Accordingly, the discovery contours for these channels in the My^-tan j3 plane deter- 
mine the boundary of the region where only the (SM-like) light CP-even Higgs boson can be 
detected at the LHC. Since the discovery contours for the LHC are at smaller values of tan (3 
compared to those accessible via the current exclusion bounds at the Tevatron, the impact 
of the tan /3-enhanced supersymmetric corrections is less pronounced in this case. We have 
studied the effect of including the dominant supersymmetric corrections, which had been 
omitted in the analyses of the production processes at the LHC, and their variation with 
the relevant parameters. Possible decays of the heavy MSSM Higgs bosons into charginos 
and neutralinos have been taken into account. We have found that the prospective discovery 
contours at the LHC are shifted by up to Atan/3 < 10. 

Based on our analysis of the sensitivities of the searches for MSSM Higgs bosons at the 
Tevatron and the LHC we have defined benchmark scenarios for the analysis of MSSM Higgs- 
boson searches at hadron colliders. They arc based on a generalization of similar benchmark 
scenarios proposed for the searches for SM-likc MSSM Higgs bosons at the Tevatron and the 
LHC. 
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